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1 Introduction 

At the present time there are known only few analytical methods enabling to solve suffi- 
ciently general sets of initial or boundary value problems for nonlinear PDEs and they are 
applied to either completely integrable PDEs or those PDEs that can be transformed 
to linear ones by means of non-local change of variables j2]. 

The majority of evolutionary PDEs used to simulate non-linear transport phenomena 
is not completely integrable. Yet it is well known that under certain conditions coher- 
ent structures formation take place during the transport phenomena occurring in open 
dissipative systems [3j and their analytical description is of great interest. 

The simplest coherent structures are represented by periodic, quasiperoidic, kink-like 
and soliton-like travelling wave (TW) solutions. In recent decades a number of effective 
methods, enabling to obtain analytical expressions for TW solutions describing coherent 
structures have been put forward [H El El 13 UHl UH 1121 UH]- But most of the papers 
dealing with this subject concentrate upon the finding out solutions, evolving in a self- 
similar mode. And only few authors look for "bi-soliton" (or a "multi-soliton") solutions 
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to those PDEs, that are not completely integrable. Being informal, we mean by this term 
solutions describing interaction of TW regimes. 

Perhaps, the most advanced study of bi-soliton solutions of PDEs that are not com- 
pletely integrable has been performed by H. Cornille and A. Gervois ^[E]. They de- 
scribed a broad class of non-linear PDEs, possessing bi-soliton solutions. Yet the method- 
ology put forward by these authors applies merely to so called "factorized" PDEs and, 
besides, further investigations revealed that their classification is not complete. T. Kawa- 
hara and M. Tanaka considered equation of the Fisher type that does not fit the clas- 
sification scheme from fH]- They succeeded in analytical description of travelling fronts, 
using the classical Hirota method [T^ . 

In this work we present the analytical description of interacting wave fronts within 
the hyperbolic generalization of Burgers equation, using the modified Hirota method. 
The structure of the study is following. In section 2 we introduce a canonical form 
of generalized Burgers equation and study its TW solutions within the classical Hirota 
method. At the beginning of section 3, using the preliminary information about the TW 
regimes, we try to obtain bi-soliton solutions using the classical Hirota method [Dodd] 
and show that it is rather impossible. So in the following part of this section we use a 
modified version of the Hirota method, enabling to achieve our goal. 

2 Hyperbolic generalization of Burgers equation and 
its TW solutions 

We consider the hyperbolic generalization of Burgers equation : 

3 

TUa + Ut + uu^ - KU^:^ = ^ai,u''. (1) 

Here u = u{t, x), r, n, are constant parameters (for physical reason r > 0, and 
K > 0). Lower indices denote partial derivatives with respect to corresponding variables. 
Assuming that the polynomial in RHS of equation has three real roots, we can rewrite 
this equation in the following equivalent form: 

TUtt + Ut + UUx + B Ux — Uxx = A M {u — S) {u — Q) , (2) 

where B, S, Q, X are constant parameters. We'll use for equation (0) the abbreviation 
GBE. 

As it was previously announced, the main goal of this study is to obtain the analytical 
description to the bi-solitons, using the Hirota method. Employment of this method, that 
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will be briefly described later on, leads to a very complicated system of non-linear algebraic 
equations that rather could not be solved directly without any additional information 
about parameters being involved into the scheme. Processing of the system of algebraic 
equations appearing within the Hirota method becomes more simple if we restrict ourself 
to finding out solutions, possessing some already known asymptotic features. We will 
assume in this study that bi-soliton solutions tend on +00 or (and) —00 to corresponding 
TW solutions. For this reason we begin with the analytical description of kink-like TW 
solutions, playing the role of asjTiiptotics. To obtain them, we use the ansatz 

n = ^. (3) 

Following the Hirota method, we put 

/ = 1 + e(p{t, x), 

where ip{t, x) is an unknown function and e is a formally small parameter. We say 
that (f{t, x) is "formally" small, since after the substitution of into (j2]) we use the 
asymptotic decomposition of the expression obtained, treating this parameter as a small 
one, but finally we put e = 1. The detailed description of this procedure can be found 
e.g. in p. 

So, we insert ansatz into (j2)) and multiply the obtained expression by This way 
a third-order polynomial with respect to e is obtained. Equating to zero the coefficient of 
e^, we get the following PDE: 

T ipttx + ^tx + B ip^x - K Lp^xx + QS\ipx = Q. (4) 

It is obvious, that solutions of the linear PDE Q can be presented in the form 

ip{t, x) = exp [ax — vt + c] , (5) 

where a, v, c are constant parameters. Inserting (0) into the ansatz Q, and equating to 
zero coefficients of e'^, /c = 1, 2, 3, we obtain the following system of algebraic equations: 

-aB + v + a^K-QS\-v^T = 0, (6) 
-V + a^{l + k) + 2Q S X + a[B + {S - Q) X] - v^t = 0, (7) 

{a + Q)X{a- S) = 0. (8) 

We see from equation (jHJ that in case when A is non-zero parameter a is equal to either 
S or —Q. 
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In the first case we get the solution of the form 

{t, x) = exp [S X — tv + ci] , 



(9) 



^^__ S[2B + S{\ + 1) + 2\Q] _ ^^^^ 

This function being inserted into Q gives rise to a kink-hke solution in case when an 
extra condition 

T[2B + S + \i2Q + S)f -2{1 + X) . . 

K = K^ = (11) 

is satisfied. As it was mentioned before, we finally put e = 1 in the formula Q. This may 
be done correctly because the parameter c in (0) is arbitrary. So we can "rescale" it as 
follows: c ^ c + log 1/e. This trick eliminates e from the formula Q. 

In the second case, i.e. when a = —Q we get the following solution of system (jUjl-lfTj): 

(P2{t, x) = exp[-Q X - tw + €2] , (12) 
-Q [2B-Qi\ + 1)-2\S] ^^g^ 

which leads to the kink-like TW solution providing that 

_ r[\i2S + Q) + Q-2Bf-2il + \) 

H — — ^ • 

We see that, generally speaking, parameters and are different. But, in contrast to 
a, V, w and q, z = 1, 2, the parameters k, characterizes equation Q and cannot be 
chosen arbitrarily. So in order that Q and (jl2p describe two different solutions of GBE, 
we should equahze and k^. It is easy to show that equalities 

_ (g + g)V(l + 3A)-8(l + A) 
lo 

take place if 
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So when k and B are given by equations (fT3|) and (fTB|) correspondingly, we get two 
different kink-like TW solutions of the same equation. 



3 Bi-soliton solutions 

3.1 Procedure based on the classical Hirota method 

We look for the solution of the following form: 

u(t, x) = F {uji, UJ2) , (17) 
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where ui = S x — vt + ci, lj2 = —Qx — wt + C2- Solution (fT7|) describes interactions of 
various TW in cases when ui and uj2 are not proportional. As it was mention before, in our 
attempts to obtain bi-soliton solutions within the Hirota method, we face the necessity 
to solve very difficult systems of non-linear algebraic equations. To make the problem 
more tractable, we should pose some conditions on the parameters to be determined. So 
assuming that solution we are looking for describes asymptotically (depending on the sort 
of interaction) one or two TW fronts, we choose parameters v and w in accordance with 
the formulae (fTTHl and |T<?|) . Of course, we can proceed this way further on and incorporate 
the formulae (fTH|l as well, but our analysis shows that in this case there is a lack of 
non-trivial solutions. 

When looking for simplest bi-soliton solution, it is instructive to choose / in the form 
of superposition of functions corresponding to the TW solutions: 



Inserting (jH)) with / given by the formula (jl8|) into equation (0) and multiplying the 
obtained expression by we get a six-order polynomial with respect to e. Equating 
to zero coefficients of e", n = 1, 2, ...6, we obtain six nonlinear algebraic equations (we 
denote them as (eql),...(eq6)), containing the parameters and products of ea;p[u;i] and 
exp[uj2]. For brevity we shall use the notation 



Calculations were carried out in several steps essentially based upon employment of the 
"Mathematica" software. Below we describe a crucial points of the procedure employed, 
enabling to repeate it. Yet it is obvious, that correctness of solutions obtained within this 
procedure can be checked by direct inspection without going into details. 

Before we proceed further on, let us note that we do not take into account solutions 
with Q = 0, S = 0, Q = —S and A = —1/3. Elementary analysis of formula Q whth 
(fTUj) and (fT5j) being taken into account shows, that in all these cases we get at most TW 
solution. 

So let us describe the procedure of obtaining bi-soliton solution based on the ansatz 



f + e {exp[uji] + exp[LU2]) + Re^ exp[uji + uj2\. 




X' 



.m 



exp[m ■ ui], = exp[n ■ 002]. 



Step 1. First of all we consider (eq6), which is the simplest one: 



QR^ S {S -Q)x^y^ X = 0. 



(19) 



From this equation we obtain that 



S = Q. 



(20) 
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Figure 1: An example of a bi-soliton solution of equation (j2Z|), described by the formula 
(EHI) for 5 = 2.5, ci = 2, C2 = -6, k = 0.5 



Step 2. Next we consider (eql), which is a linear one with respect to variables x, y. Taking 
into account and equating to zero coefficients of x and y, we obtain the following 
system: 



-2(l + A)-4K + r [4 52 + 45 5(1 + 3A) + 52(1 + 6A + 9A2)J =0, (21) 
2(l + A)+4K + r \-4B^ + 4B S{1 + 3X) - S^1 + 6X + 9X^)] =0. (22) 



Summing up these two equations we get the equation 



8BSt (1 + 3A) = 0. 



(23) 



To fulfill (j2ni), we put B = 0. Returning then to equation (PT|) we obtain the 
following expression for k: 

-2(1 + A) +r52(l + 6A + 9A2) 



K 
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Step 3. With these conditions (eql), (eq2) and (eq6) nullify while subsystem (eq3)-(eq5) 
becomes as follows: 



S^xy{x-y){l + 3\) 2S\{1- R){l + 3\)+R 



RS'-^ xy (x 



3 A) 



0, 
0, 



R' S'x' y' (x-2/)(l + 3A) = 0. 



(24) 
(25) 
(26) 



The only non-trivial solution to ()24j ) - (PBj) corresponding to bi-solitons is i? = r = 0. 
It can be easily checked starting from equation (j^^ that we do not obtain any other 
solution choosing solution r = instead of 5 = in the second step. 
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Thus, we became convinced that there is impossible to obtain a bi-sohton solution 
satisfying (j21), using the classical Hirota method. Yet as a by-product we get the following 
result: 

Theorem 1. The Burgers equation 

Ut + uu^ — K = — (1 + 2 n {u^ — S*^^ (27) 
possess a bi-soliton solution 

exp{Sx -vt + ci) -exp{-Sx -vt + C2) 

" ^' = S ^ 7^ ^ 7 (28) 

1 + exp [b X — V t + ci) + exp {—b x — vt + C2) 

where v = —S"^ (1 + 3 k) . 

Example of bi-soliton solution describing interaction of two wave fronts is shown in 
Fig. □ 

3.2 Procedure based on the modified Hirota method 

One can succeed in obtaining bi-soliton solutions to ^ by slightly modifying Hirota 
method. Let us consider the following ansatz: 

uit, x) = ^, (29) 

in which we put 

g = e[a exp (cui) + P exp (0^2)] + e'^ A exp {ui + UJ2) , (30) 



f = 1 + e [exp (ui) + exp {UJ2)] + Rexp {ui + UJ2) ■ (31) 

So we insert ()29|) into multiply the resulting equation by and obtain this way a 
six-order polynomial with respect to e. Equating to zero coefficients of e'', k = 1,2, ...6, 
we get a system of non-linear algebraic equations which we again denote by (egl)-(eg6). 
Using the "Mathematica" software we succeeded in obtaining a non-trivial solustion of 
these equations. A sketch of the procedure employed is presented below. 

Step 1. We consider (eq6) which reads as follows: 

-A (A + QR) {A- S R)x^y'^ = 0. 

In order to satisfy this equation, we put A = RS. 
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Step 2. Next we consider (eql). Nullifying the coefficient of we obtain the equation 



Q^(3 {-2{l + 2K) + r{AB^-ABQ + Q^) + {Q + 2SfX\+ 
+X [-2 (l -QV) + AQ S t - a B t {Q + 2 S)]} = 0. 

This equation will be satisfied if we choose (3 = 0. With this choice we consider 
(eq5) and, nullifying coefficient of obtain: 

{S -a)UQS -AS^X + tQ^ (1 + A)^ + 4gV {l + X){SX-B) + 



+2Q'^{1-2k + X + 2B'^t-ABSXt + 2 S^X'^t 



0. 



This equation will be satisfied if S* = a. Equating to zero the remaining term, we 
get: 



jR^S^x^y^ {-2 - Ak + AB\ + AB St + S\ + {2Q + Sf X' 



(32) 



-2 A 



-1 + 2Q S T + S\ + 2 B {2Q + S) t]] = 



Equation ()32|) is satisfied if n is expressed by the formula 
Step 3. Next we consider (eq3). The coefficient of x'^y reads as follows: 

-^g(i?-l)^(Q + ^)'(l + 3A)r(4 5 + 2^-g + Ag + 2A5) 



In order to nullify it we put 



B 



2S-Q + X{Q + 2S) 



(33) 



The remaining coefficient of (eq3) to be nullified is that of xy"^. Taking (|33|) into 
account we can write it as follows: 

-QiR-l)SiQ + S + QX + 2SX). 

Analysis of formula shows that option R = 1 leads to a one-soliton solution so 
we nullify this expression assuming that 

Q + S 



X 



Q + 2S' 



(34) 



Step 4. Taking into account all the above conditions and repeating the procedure we con- 
clude that only (eq2) and (eq4) are nonzero: (eq2) is given by the formula 

{R-l)S^xy[2Q + AS + AQ^T + 8Q^ST + 5QS^T + TS^] 

2{Q + 2Sf ' 



Figure 2: Creation of two wave fronts "from nothing", described by the formula (jHfjj) for 
S = 1,Q = -1.9, ci = 10, C2 = -6, R = 




Figure 3: Interaction of locahzed TW solution and kink described by the formula ()36j) for 

5 = 1, g = 2, ci = 2, C2 = -3, i? = 1500 

while (eq4) is as follows: 

RQ^R-1) 5 V y"^ [2Q + 4S + 4Q^T + 8Q^S T + 5Q S^r + r S^] 



Both of them are nullified if 



2{Q + 2Sf 



2{Q + 2S) 
(Q + S) {2Q + Sf' 



(35) 



And this way we obtained the bi-soliton solution to GBE. Elementary but slightly 
tiresome analysis shows that system (GBE) does not admit any other bi-soliton solution 
described by the formula ()29|) . The only exception is the solution connected with the 
already obtained solution by means of the transformation S (-^ Q, Q ^ —S, a = 
0, /? = —Q. But existence of the second solution is a mere consequence of the arbitrariness 
of the choice of what is Ui and what is UJ2 in the formulae ()30 p - (j31|) . It is obvious then, 
that this way we don't get any other independent solution. 
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Let us summarize the result obtained as the following statement. 



Theorem 2. // 



r 



2{Q + 2S) 



and 



Q + S 

'Q + 2S' ■ {Q + S){2Q + Sy 
Q + 2S 



8 (Q + S) 

then equation ^ admits the following bi-soliton solution 

S exp (cji) [1 + Rexp {UJ2)] 



u{t, X 



where 



ui = S 



X 



1 + exp (ui) + exp {U2) + Rexp {ui + U2 
t Q{2Q + Sy 



4 Q + 2S 
S, Q, R, Ci, C2 are arbitrary constants. 



Ci, LJ2 = -Q 



x--{2Q + S) 



C2, 



(36) 



Patterns of the interacting wave fronts described by the formula ()36p are shown in 
Figs. H El 



4 Concluding Remarks 

Thus, we have obtained bi-soliton solutions for classical Burgers equation with cubic 
non-linearity and for the GBE. To our knowledge, such solutions for the equation Q 
have been obtained for the first time. In this study we looked for the simplest bi-soliton 
solutions described by the formula (j^^ . It is obvious that rational solutions of this sort 
can be unlimitedly generalized. However it should be noted, that one faces very non- 
trivial technical problems when incorporating additional terms into (j3(J|) - (|31|) . And for 
successful application of the above method it is desired to have some extra information 
about the possible values of the parameters. Beside the asymptotic analysis, one can 
incorporate another methods, such as the Painleve analysis jlH] and study of conserved 
quantities [TU] . 

Let us mention in conclusion that solutions obtained in this work describe merely inter- 
action of kinks. It would be interesting to obtain an analytical description of interacting 
TW waves of another sort. In paper [20] there were obtained periodic, soliton-like and 
some other solutions to GBE. All of them have been obtained within the methods different 
from that used in this paper. Consequently a successful investigations of the interactions 
of TW of different types should be based on employment of more sophisticated ansatzes. 



10 



References 



[1] R.K. Dodd, J.C Eilbek, J.D. Gibbon, H.C. Morris, Solitons and nonlinear wave 
equations, Academic Press, London, 1984. 

[2] G. Bluman, S. Kumei, Symmetries and Differential Equations, Springer, New York, 
1989. 

[3] P. Glansdorf, I. Prigogine, Thermodynamics of structure, stability and fluctuations, 
Wiley, New York, 1971. 

[4] W. Fushchich, I.Tsyfra, On a Reduction and Solutions of Nonlinear Wave Equations 
with Broken Symmetry, J. Phys. A: Math. Gen., 20 (1985), L45-L48. 

[5] W. Fushchych, W. Shtelen, N. Serov, Symmetry Analysis and Exact Solutions of 
Equations of Nonlinear Mathematical Physics, Kluver Publishers, Dordrecht, 1993. 

[6] P. Clarkson, M. Kruskal, New Similarity Reductions of the Boussinesq Equation, J. 
Math. Mech., 18 (1989), 1025-1042. 

[7] P. Clarkson and L. Mansfield, Symmetry Reductions and Exact Solutions of a Class 
of Nonlinear Heat Equations, Physica D 70 (1993), 250-288. 

[8] P. Olver and P. Rosenau, The Construction of Special Solutions to Partial Differential 
Equations, Phys. Lett., A 114 (1986), 107-112. 

[9] P. Olver, E. Vorob'jev, Nonclassical and Conditional Symmetries, in: CRC Handbook 
of Lie Croup Analysis of Differential Equations, Volume HI, ed. by Nail Ibragimov. 
- CRC Press Inc, 2000. 

[10] E.Fan, Multiple Travelling Wave Solutions of Nonlinear Evolution Equations Using a 
Unified Algebraic Method, Journ of Physics A: Math and Gen 35, (2002), 6853-6872. 

[11] A.Barannyk, I.Yurik, Construction of Exact Solutions of Diffusion Equation, Pro- 
ceedings of Institute of Mathematics of NAS of Ukraine 50, Part I (2004), 29-33. 

[12] A. Nikitin, T. Barannyk, Solitary Waves and Other Solutions for Nonlinear Heat 



Equations, preprint, http://arxiv.org /abs/math-ph/0303004[ 2003. 



[13] V.A. Vladimirov and E. V. Kutafina, Exact Travelling Wave Solutions of Some Non- 
linear Evolutionary Equations, Rep. Math. Physics 54 (2004), 261-271. 



11 



[14] H. Cornille, A. Gervois, Bi-soliton Solutions of a Weakly Nonlinear Evolutionary 
PDEs with Quadratic Nonlinearity, Physica D, 6 (1982), 1-50. 

[15] H. Cornille, A. Gervois, Connection Between the Existence of Bisolitons for Quadratic 
Nonlinearities and a Factorization of the Associate Linear Operator, J. Math. Phys., 
24 (1983), 2042-2055. 

[16] T. Kawahara, M. Tanaka, Interactions of Travelling Fronts: an Exact Solutions of a 
Nonlinear Diffusion Equation, Phys. Letters, 97A (1983), 311-314. 

[17] A.S. Makarenko, Mathematical modelling of memory effects influence on fast hydro- 
dynamic and heat conduction processes, Control and Cybernetics 25 (1996), 621-630. 

[18] M. Tabor, Chaos and Integrability in Nonlinear Dynamics: An Introduction, Wiley, 
New York, 1989. 

[19] G.Bluman, Connections Between Symmetries and Conservation Laws, Symmetry, 
Integrability and Geometry: Methods and Applications, 1 (2005), Paper Oil, 16 p 
( |http:/ /www.emis.de/journals/SIGMA/2005/Paper011/index.htmi] ). 

[20] V.A. Vladimirov and E.V. Kutafina, Analytical Description of Coherent 
Structures within the Hyperbolic Ceneralization of Burgers Equation, preprint 
http: / /arxiv.org/abs/nlin .PS/0 607013t 2006. 



12 



